Direct measurement and reconstruction of nonclassical features of twin beams 
generated in spontaneous parametric down-conversion 
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Correlations in twin beams composed of many photon pairs are studied using an intensified CCD 
camera. Joint signal-idler photon-number distribution and quantum phase-space quasi-distributions 
determined from experimental data have nonclassical features. 
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Quantum mechanics interprets nonlinear optical pro- 
cesses as physical effects composed of many elementary 
'quantum' events in which one (several) photon is anni- 
hilated and several (one) photon emerge Ml. This behav- 
ior of nonlinearly interacting physical systems is com- 
pletely unusual in classical physics and lies in heart of 
nonclassical-light generation. Generation of correlated 
photon pairs in the process of spontaneous parametric 
frequency down-conversion is probably the most com- 
mon example 0- Fundamental experiments performed 
for the first time with photon pairs have shown that these 
nonlinear elementary 'quantum' events have even an 'in- 
ternal structure'. Measurement using photon pairs in 
Hong-Ou-Mandel interferometer has revealed this struc- 
ture demonstrating that photons in one photon pair (pho- 
ton twins) are generated within a sharp time window 
typically of several tens of fs . Later experiments have 
even shown that entangled photon pairs are the same 
fundamental entities as single photons. Similarly as a 
single photon can interfere only with itself, an entan- 
gled photon pair interferes only with itself Q. These 
completely unusual properties of photon pairs have been 
experimentally tested in numerous experiments with the 
same qualitative conclusion - fundamental laws of quan- 
tum mechanics are a solid basis for the explanation of 
obtained experimental results. Among others, experi- 
mental confirmation of violation of Bell inequalities us- 
ing photon pairs has excluded neoclassical theories with 
local hidden variables as a right tool for the description 
of Nature. Photon pairs are also an indispensable tool 
in quantum teleportation quantum cryptography or 
dense coding. 

The use of intense femtosecond pump fields together 
with availability of new materials with larger nonlincar- 
ities have opened a new area in investigation of twin 
beams. Nowadays even beams containing many photon 
pairs generated in a sharp time window can be obtained. 
The first experiments confirm in agreement with quan- 
tum theory that photon pairs inside such beams behave 
as independent entities. Coherence properties of twin 
beams reflect those of the pump beam. Moreover, it has 
been shown that photonpairs can have their origin also 
in stimulated emission [6( (laser-like generation of twin 
beams is sometimes mentioned). 

These fundamental properties of photon pairs then 
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FIG. 1: Scheme of the experiment. The inset shows accumu- 
lated picture after 240,000 frames. 



determine statistical properties of 'more intense' twin 
beams. In this letter, we report on experimental de- 
termination of photocount statistics of twin beams 
using an intensified CCD camera (iCCD). An alterna- 
tive approach to determine photocount statistics is based 
on homodyne detection and has revealed correlations in 
photon numbers of two fields comprising a two-mode 
squeezed state 0- In our case, even raw experimental 
data show nonclassical properties of twin beams. The re- 
constructed photon-number distribution then enables to 
determine (radially symmetric) joint phase-space quasi- 
distributions of the twin beams. 

In the reported experiment sketched in Fig. ^ we 
use the second-harmonic of an amplified femtosecond 
Ti: sapphire laser system as pump source for the gen- 
eration of down-converted light in a 5-mm long Li:I03 
nonlinear crystal. Laser system operates at 800 nm and 
yields a train of 200 fs pulses with a repetition rate con- 
trolled by the amplifier and set to 11 kHz in our case. 
Their second-harmonic is produced in a 2-mm long BBO 
crystal and pulses of energy up to 1.3 /jJ are obtained. 
Down-converted photons emerge from a crystal at a cone 
layer with vertex half-angle of 31 deg (see Fig. QJ. One 
part (signal) of the cone layer is captured by the photo- 
cathode of the camera placed 10 cm behind the crystal. 
The corresponding opposite part of the cone layer (idler) 
is led to the camera after reflection on a high-reflectivity 
mirror placed as close to the crystal as possible to mini- 
mize path difference of the two beams. Prior to entering 
the camera, both beams are filtered by two edge filters 
(high-pass above 750 nm) and a 20 nm wide (FWHM) 
bandpass filter centered at 800 nm. The filters reduce the 
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noise coming from fluorescence in the nonlinear crystal 
and stray light in the laboratory to an acceptable level 
and bandpass filter defines the thickness of the cone layer. 
In the software of the camera, three regions of interest 
are defined, two for the signal and idler strips and a third 
one that serves for monitoring the noise level (see inset 
in Fig. The overall quantum efficiency of the iCCD 
camera has been found to be 9.5%. 

Quantum state of twin beams at the output plane of 
the crystal can be described by the following statistical 
operator psi 0: 
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\ n s)s (\ n i)i) denotes Fock state with ns (ni) photons 
and then p(ns,n{) means the joint signal-idler photon- 
number distribution. The down-converted photon pairs 
then undergo several loss mechanisms before they are reg- 
istered. Taking into account losses in the signal and idler 
paths (expressed using effective transmittances Tg, Tj), 
quantum efficiency of the camera (77), and noise due to 
other light sources and internal noises of the camera (D), 
the measured probabilities f(cs,cj) of having cs detec- 
tions in the signal strip and cj detections in the idler strip 
at the camera are determined as: 
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This formula holds provided that the number of photons 
detected by the camera is much lower than the number 
of active pixels 

An example of the measured photon-number distribu- 
tion f(cs, cj) is given in Fig.[2t- A correlated character of 
twin beams is clearly visible in the graph in Fig. [SJd show- 
ing the difference between the measured photon-number 
distribution /(cs, ci) and that one composed of two inde- 
pendent Poissonian distributions with mean values equal 
to those of the measured signal and idler fields. We can 
see that elements lying on diagonal or near diagonal are 
enhanced whereas those lying far from diagonal are sig- 
nificantly suppressed. This is a direct experimental man- 
ifestation of the fact that signal and idler photons are 
generated in pairs. 

Covariance C p of the signal and idler photon numbers 
determined along the formula [f0| 
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FIG. 2: (a) Measured joint signal-idler photon-number distri- 
bution / as a function of the signal (ns) and idler (nj) photon 
numbers, (b) difference between the measured photon-number 
distribution / and the distribution given by direct product of 
two independent Poissonian distributions with means equal to 
those of experimental data. Solid line denotes zero contour. 



equals 0.0435 ±0.008 for the data shown in Fig.0 Type 
of photon-number statistics can be judged according to 
the value of coefficient K, 



K 



(n) 2 (n) 



(5) 



K = 1 for Poissonian statistics, K = 2 characterizes 
Gaussian statistics, and K < 1 for a Fock state. The 
experimentally determined marginal signal-field (idler- 
field) photon- number distribution has Ks (Ki) equal to 
0.997 ± 0.030 (0.994 ± 0.030), i.e. the statististics can 
be considered to be Poissonian within the experimental 
error. Statistics of photon pairs deviate from the Poisso- 
nian ones if the number of independent modes constitut- 
ing the field is small . 

Taking into account the fact that the measured 
marginal signal and idler photon-number distributions 
are Poissonian, the measured nonzero value of covariance 
C p also indicates nonclassical properties of detected twin 
beams. We note that a classical coherent field (with Pois- 
sonian statistics) cannot have any correlations in photon 
numbers at two distinct spatial points. 

Pair character of twin beams is clearly revealed in 
the original joint signal-idler photon-number distribution 
p(ns,m) characterizing the light field as it occurs at 
the output plane of the crystal. The method of maxi- 
mum likelihood estimation has proven to be extraordi- 



3 



a) 




110 ~120 
120 



b) 




50 60 70 



100 110 120 



FIG. 3: (a) 3D- and (b) topograph of reconstructed joint 
signal-idler photon-number distribution p'°°^ as a function of 
the signal (ns) and idler (nj) photon numbers; area in graph 
(b) bordered by the solid line contains probabilities violating 
classical inequality in Eq. r)T s = 0.0539, rfTi = 0.0415, 
D = 0.75. 



narily suitable in finding the inversion to the relation in 
Eq. [3 The joint signal-idler photon-number distribution 
P^°°H n s, n i) that minimizes the Kullback-Leibler diver- 
gence from experimental data can be found using the 
iterative Expectation-Maximization algorithm [T^: 
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The symbol p^ n \ns,ni) denotes the joint signal-idler 
photon-number distribution after the rt-th step of iter- 
ation, p(°\ns,ni) is an arbitrary initial photon-number 
distribution. 

The joint signal-idler photon-number distribution 
p(°°'(ns,ni) at the output plane of the crystal deter- 
mined using Eq.|^l shows strong correlations between the 
signal photon numbers and idler photon numbers (see 
Fig. |3J. Values of the parameters nTs, rjTj, and D nec- 
essary for reconstruction have been determined from an 
independent measurement. Covariance of the signal and 
idler photon numbers equals C p « 0.98 after the re- 
construction and shows the ability of the reconstruction 
method to reveal the original pairing of photons. This 



reconstruction is not perfect due to the impossibility to 
describe precisely all noises occurring in the experiment. 
Mixing the signal and idler fields together, the resul- 
tant photon-number distribution preferring even photon- 
numbers also reflects these correlations, as has been ob- 
served in |l3j | using a special photon-number resolving 
detector. 

Any photon-number distribution pins, nj) originating 
in a classical field has to fulfil the following inequality 
0: 

p(ns, ni) < -$-r exp(-ns) -^7 exp(-nj); (7) 

i.e. this inequality represents a criterion of nonclassical- 
ity. The reconstructed probabilities lying inside the bold 
curve in Fig. 3b exhibit violation of this inequality. The 
violation is a consequence of the generation of photons in 
pairs (probabilities tend to be concentrated towards the 
diagonal where they reach greater values). 

Experimental marginal signal and idler photon-number 
distributions derived from the reconstructed joint signal- 
idler photon-number distribution are Poissonian (Kg — 
0.997, Kj = 1.000). This means that also the probabil- 
ity distribution of generated photon pairs is Poissonian. 
Poissonian statistics reflects the fact that photon pairs 
can be generated into many independent modes distin- 
guishable in space and time. In case of our experimen- 
tally obtained twin beams having typically tens or hun- 
dreds of photon pairs, each pair occupies its own mode 
with a high probability. 

The reconstructed joint signal-idler photon-number 
distribution p*- 00 ' enables to determine the joint distri- 
bution of integrated intensities Ws and Wj of the signal 
and idler fields at the output plane of the crystal; 
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Inverting photodetection equation we have for the 
joint signal-idler integrated-intensity distribution 
P(Ws, Wi,s, M) related to s~ orderin g of field operators 
and assuming M independent modes f|l5|. chap. 4): 
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— 1 are Laguerre polynomials and T means the gamma 
function. The parameter s equals -1, 0, and 1 for anti- 
normal, symmetric, and normal ordering of field opera- 
tors, respectively. 

Since both signal and idler fields are phase inde- 
pendent (they are generated in a multimode sponta- 
neous process), the joint signal-idler phase-space quasi- 
distribution $ (ag , a/ , s, M) = <&( \ots | , \ai \ , s, M) related 
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FIG. 4: Joint signal-idler quasidistribution 4>(|qs|, \ai\, s — 
0.3, M — 50) determined from (a) experimental data and 
(b) assuming an ideal signal-idler photon-number distribution 
p(ns,ni) = (Sns.nj exp(-/i)/i" s /n s !, fi = 14.5. 



to s— ordering can be simply determined: 

§(\a s \,\<xi\,s,M) = \p{\a s \ 2 ,\ai\*,s,M) (10) 



using integrated-intensity distribution P. A nonclas- 
sical field manifests itself by negative values in quasi- 
distributions with s > sq , where the value of so is a char- 
acteristic of the state. Quasi-distribution with s = 0.3 
obtained from a reconstructed photon-number distribu- 
tion p°°{ns,ni) and plotted in Fig. 0Ji has a typical 
structure of correlated states composed of regions with 
positive and negative values. For comparison, the same 
quasi-distribution $ideai (| ens |, \a/\,s,M) belonging to an 
ideal (perfectly correlated) Poissonian joint signal-idler 
photon- number distribution is shown in Fig.QJx Regions 
with neg ative values are the evidence of nonclassical na- 
ture [16| of twin beams composed of many photon pairs. 

The iCCD camera is also able to reveal spatial cor- 
relations of photon pairs |17| that can be characterized 
by entanglement area of two photons in a photon pair. 
Spatial correlations can then be used for further studies 
of photon-number statistics. These investigations are in 
progress. 

In conclusion, we have experimentally demonstrated 
that twin beams composed of many photon pairs are 
of nonclassical origin that is manifested both in joint 
photon-number distributions and joint quantum phase- 
space quasi-distributions. 
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